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Abstract 

We present derivations of the expressions for the spherical volume averages of static electric and 
magnetic fields that are virtually identical. These derivations utilize the Coulomb and Biot-Savart 
laws, and make no use of vector calculus identities or potentials. 
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I. INTRODUCTION 



The average of static electric or magnetic fields over a sphere has been used to obtain 
important results, such as the macroscopic electric field inside a dielectrioi and the presence 
of 5-function electric (magnetic) fields at the position of electric (magnetic) point dipoles.- 
In both electric and magnetic cases, if the sources of the fields are outside the averaging 
sphere, the average field equals the value of the field at the center of the sphere, whereas for 
sources inside the averaging sphere, the average electric and magnetic fields are proportional 
to the electric and magnetic dipole moments of the sources, respectively. 

Textbooks typically only treat the electric field case.- Only a handful of authora^i^ treat 
both electric and magnetic cases, probably because the magnetic field case is considered to be 
"tough" by undergraduate standards,- since the derivations typically employ the magnetic 
vector potential and vector calculus identities. This paper presents derivations of both the 
electric and magnetic field cases that are virtually identical and are elementary, in the sense 
that they rely on Coulomb's and the Biot-Savart laws and make no use of vector calculus 
identities or potentials. 

Coulomb's law states that the electric field for a unit point charge is 

G{r) = fc^, (1) 

where r is a displacement vector from the point charge, r = |r | and k = (47reo)^^ and 1 in 
SI and Gaussian units, respectively. In terms of G, the electric field the electric field E at 
point s due to a static charge distribution p(r') is 

E{s) = j dr' p{r') G{s - r') = - j dr' p{r') G{r' - s). (2) 

Define Gr{r') to be the average ofG over a spherical surface of radius r around a point r', 
as shown in Fig. 1. An apparently little-known result that is critical in following discussion 
is (see the appendix for derivations) 




where B is the Heavyside step function.^ Eq. ([3]) implies that the average of the electric 
field over the surface of the sphere for a point charge outside the sphere (r' > r) equals the 
electric field at the center of the sphere, whereas a point charge inside a sphere (r' < r) 
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contributes zero to the average of the electric field over the surface of the sphere. In a sense, 
this result is the antithesis of the integral form of Gauss' law/ but it is not as general as 
Gauss' law because it applies only to spherical surfaces. 

II. VOLUME AVERAGES OF STATIC ELECTRIC FIELDS 

Let {E)ii to be average of the electric field E about a spherical volume of radius R and 
Er{Q) be the average electric field over a spherical shell of radius r, both centered around 
the origin. Then, 

For a given charge distribution p{r'), Er{0) can be obtained by setting s = and averaging 
both sides of Eq. ([2]) over a spherical shell of radius r, yielding 

E;{Q) = - j dr' p{r')G;{r'). (5) 

A. Single point charge 

We first examine the case of a single point charge. Because electric fields obey the 
principle of linear superposition, what holds for a single point charge is generalizable to 
cases of many charges and/or continuous charge distributions. Substituting the charge 
distribution for a point charge q at d, p{r') = q6{r'—d), into Eq. ([5]) gives Er{0) = —qGr{d). 
Substituting this into Eq. (jlj) and using Eq. ([3]) gives 

{E)n = -^^iyrr'e{d-r). (6) 

If the point charge is outside the averaging sphere of radius R, then d > R and the integral 
in Eq. (El) equals R^/3. Hence, (E^)^*"*^ = —kqd/d^ = E{0), the electric field at the center 
of the sphere due to the point charge. On the other hand, if the point charge is inside 
the averaging sphere, then d < R and the integral in Eq. ([6]) equals d^/3, in which case 
(E)^^^ = —kp/R^, where p = qd is the dipole moment of the point charge relative to the 
center of the sphere. 
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B. Arbitrary charge distribution 



We now obtain these results more rigorously for an arbitrary charge distribution p{r'). 
Substituting Eqs. P and © into Eq. (HD yields, 

-n- JO \_Jr'=r ' 

where the effect of the G function in Eq. ([3]) is to restrict the r' integration to the region 
\r'\ > r. We now consider separately the contribution of the charges outside and inside the 
sphere. 



(7) 



1. Sources outside the sphere 

The contribution due to charges outside the sphere of radius R corresponds to dr' 
in Eq. ([71), which makes the term in square parentheses independent of r. The integration 
over r gives R^/3, yielding 

= - r dr' = ^;(-t)(o) (8) 

Jr'=R 

the electric field at the origin. Therefore, in general, (-E)^*^*^ = E^"^^^ at the center of the 
sphere. 



2. Sources inside the sphere 

Sources inside the sphere of radius R correspond to \r'\ < R 'm. Eq. (|7]). This gives 

/^v(in) 3 /-^ , 2 /-^ , , kpir') r' 



D3 / r'^ 
-n- -/r'=0 ' Jo 



where p = J dr' p{r') r' is the dipole moment relative to the center of the sphere, and we 
have used dr jf dr' = dr' dr, as shown in Fig. 2. 
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III. MAGNETIC FIELD CASE 



A static magnetic field B is related to the charge current density J by the Biot-Savart 
law, 

B{s) = J dr' J{r') X G{s - r'), (10) 

where k = /io/(47r) or c"~^ for SI or gaussian, respectively. Comparison of Eq. fllOj) with 
Eq. ([2]) shows that the derivation for the magnetic field case can be copied wholesale from 
the electric field case by substituting E ^ B and p ^ Jx at every step. (The only point of 
caution is that the cross product anti-commutes, so care must be taken to preserve the order 
of J and G or r'.) Using these substitutions in Eq. (IHl) gives, for current sources outside 
the averaging sphere, 

= - r dr' ^-^^^^ = BiO), (11) 

Jr'=R ^ 

the magnetic field at the center of the sphere. For current sources inside the averaging 
sphere, using the substitutions in Eq. ([9]) gives 

where m = | J dr' r' x J{r') is the magnetic dipole moment.- 



APPENDIX A: DERIVATION OF EQ. ([3]) 

By definition, 

G'r{r') = ^ I dA' G{r"), (Al) 

where the integral is over the spherical shell of radius r around r'. Choose r' = r'z, 
where z is the unit vector in the 2;-direction. By azimuthal symmetry, Gr{r'z) = Gr^z^ = 
(47rr^)"^ J dA" Gz{r")z, where Gz is the 2;-component of G. By dividing the spherical surface 
into strips dA" = r'^sinOdO (see Fig. 3), Gr^z = | Iq dO sin 6* Gz{r"). Letting /i = cos 6*, and 
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using Gz{r'' 



kr'!,/\r"\^ gives 



Gr 



(j.2 _|_ _|_ 2r'rn)^/'^ 




rfi + r' 



J.2 _j_ _j_ 2r'rfi)^/'^ 



r'2 + 2r'r/i 



for r' < r; 
for r' > r, 



k d 
'2d? 



I I I I ^ I 



(A2) 



which imphes Eq. (jS]). 

This resuh can also be obtained by using the well known property that the average of 
an electrostatic potential of a point charge Q over the surface of a sphere of radius r equals 
the potential at the center if the charge is outside the sphere, and (5/(47reor) if the charge 
is inside.- Stated mathematically, if V{r") = k/\r"\, then the average of V over a spherical 
shell of radius r around r' is 



VJr') = k 



for Ir'l > r; 



for Ir'l < r. 



Averaging over spherical shells on both sides of the relation G{r') 
G~.{r') = -VVr{r'). Using Eq. in this gives Eq. (El). 



(A3) 

-V'V{r') yields 
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^ The Heavyside function is defined to be ©(x) = 1 for x > 0, 0(x) = for x < and 0(0) = ^. 
Note the difference between the average of a vector field G over a surface A, dA G/\A\ (where 
1^1 is the area of surface „4), and the flux of G through the surface, J^^dA h ■ G (where n is a 
unit vector perpendicular to area element dA). This paper deals with the former, while Gauss' 
law pertains to the latter. 
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Figures 




FIG. 1: Definition of Gr{r'): the average of the electric field of a unit charge about a spherical 
surface of radius r, centered at r' relative to the charge. 
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FIG. 2: Integration region for the limits dr dr' (grey region). Integration over this region 
can also be written as dr' Jq dr. 
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FIG. 3: Spherical surface over which Gz is averaged to obtain z- The hatched strip indicates 
the area element dA" on the sphere with constant Gz{9) = kr'!,/\r"\^ where r'^ = r' + rcosO and 



V + + 2rr' cos ( 
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